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Abstract 

We consider the motion of a planar rigid body in a potential flow with circulation and subject 
to a certain nonholonomic constraint. This model is related to the design of underwater vehicles. 

The equations of motion admit a reduction to a 2-dimensional nonlinear system, which is in- 
tegrated explicitly. We show that the reduced system comprises both asymptotic and periodic 
dynamics separated by a critical value of the energy, and give a complete classification of types of 
the motion. Then we describe the whole variety of the trajectories of the body on the plane. 
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1 Introduction and outline 



In this paper, we consider the motion of a planar rigid body surrounded by an irrotational perfect 
fluid. It ia assumed that there is a given amount n of circulation around the body, and that the body is 
subject to a certain nonholonomic constraint, which models a very effective keel or fin attached to the 
body. In the absence of circulation this system was termed the hydrodynamic Chaplygin sleigh in [IT], 
since in the absence of the fluid, the nonholonomic constraint models the effect of a sharp blade in the 
classical Chaplygin sleigh problem [7J which prevents the sleigh from moving in the lateral direction. 

The hydrodynamic Chaplygin sleigh in the presence of circulation was recently considered in [12], 
where it was shown that it is an LL system on a certain central extension of SE{2) by R 3 , where the 
cocycle encodes the effects of the circulation on the body. 

Our model for the nonholonomic constraint, which respects the Lagrange-D'Alembert principle, 
has been considered in the aerospace engineering community [23], while robotic models for underwater 
vehicles taking into account the effects of circulation were considered in [15J. 

History of the Kirchhoff equations. The motion of a rigid body in a potential fluid in the absence 
of external forces was first described by Kirchhoff [16]. His crucial observation was that the effect of 
the fluid on the body could be described entirely in terms of the added mass and added inertia terms, 
which depend on the geometry of the body only, and can be calculated analytically for a wide class 
of body shapes. Kirchhoff 's solution was extended to the case of rigid bodies moving in potential flow 
with circulation by, among others, Chaplygin [6] and [17], who derived the equations of motion for this 
system, provided an explicit integration in terms of elliptic functions, and described qualitative features 
of the dynamics, such as periodic motions. In recent years, these ideas have been extended to the case 
of rigid bodies interacting with point vortices [261 [2], vortex rings [27] and other vortical structures, 
and they have been used to describe underwater vehicles [18J and the motion of bio-organisms [U [5]. 
A comprehensive overview of the history of these equations can be found in [3]. 

Contributions of this paper. We show that the hydrodynamic Chaplygin sleigh with circulation 
is a new example of a completely integrable nonholonomic system: we discuss qualitative features of 
the dynamics, and we explicitly integrate the reduced equations of motion. 

Contents of the paper. In section [2] we recall Kirchhoff 's equations for a planar rigid body moving 
in a potential fluid and consider the Chaplygin-Lamb equations for the motion of the body in the 
presence of circulation. For the purpose of completeness, the added masses are computed explicitly for 
a body of elliptical shape. 

In section [3j the reduced equations of motion on the coalgebra se(2)* are written explicitly in 
terms of the added masses and the coefficient a that depends on the position and orientation of the 
fin (or blade) on the body. Again, for the purpose of completeness, we give explicit formulas for the 
parameters that enter the equations for a body having an elliptical shape. It is shown that the presence 
of circulation around the body adds new features to the motion of the hydrodynamic Chaplygin sleigh, 
and that there exists a critical value of the kinetic energy of the fluid-body system that divides periodic 
from asymptotic behavior. Indeed, for small, subcritical energy values, the reduced dynamics are 
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periodic. In this case the circulation effects drive the dynamics (via the so-called Kutta-Zhukowski 
force). On the contrary, if the initial energy is supercritical, then the inertia of the body overcomes 
the circulation effects and it evolves asymptotically from one circular motion to another, where the 
limit circumferences have different radius. This resembles the motion of the body in the absence of 
circulation treated in [11]. Moreover, we identify 7 regions in the reduced phase space that yield distinct 
qualitative dynamics of the motion of the sleigh on the plane. 

In section [4] the reduced equations of motion for the hydrodynamic Chaplygin sleigh with circulation 
are integrated explicitly for a generic sleigh. The form of the solution varies according to the energy 
regime. If the energy is subcritical, the solution is a quotient of trigonometric functions and we give 
an explicit expression for the period in terms of the energy. We also give a closed expression for the 
angular part of the monodromy matrix that is involved in the reconstruction process. This formula 
allows us to show that the qualitative behavior of the sleigh on the plane is very sensitive to initial 
conditions for energy values that are slightly subcritical. 

On the other hand, for the critical value of the energy, the solution of the reduced equations is 
given as a rational function of time. The closed expression for the solution is used to show that 
the motion of the sleigh on the plane in this case is bounded and evolves from one circular motion 
to another, where the limit circumferences have equal radius. Finally, the solution of the reduced 
equations for supercritical energies is given as a quotient of hyperbolic functions. In this case we also 
express the distance between the centers of the limit circumferences in terms of integrals generalizing 
the Beta-function. 

In Conclusions we motivate a further study of the hydrodynamic Chaplygin sleigh in the presence 
of point vortices. 

2 Preliminaries: Fluid-structure interactions 

In this section, we give an overview of the Kirchhoff equations describing the dynamics of a rigid body 
in potential flow, and of the Chaplygin-Lamb equations dealing with rigid bodies in the presence of 
circulation. Most of the material covered in this section can be found, for instance, in [13] as well as 
in the classical works of Lamb [T7] and Milne-Thomson [20] . 

2.1 Kinematics of rigid bodies and ideal fluids 

Following Euler's approach, consider an orthonormal body frame {Ei, E2} that is attached to the body. 
This frame is related to a fixed space frame {ei,e2} by a rotation by an angle 9 that specifies the 
orientation of the two dimensional body at each time. We will denote by x = (x, y) E M 2 the spatial 
coordinates of the origin of the body frame (see Figure [I]). The configuration of the body at any time 
is completely determined by the element g of the two dimensional Euclidean group SE{2) given by 



We will often denote the above element in g E SE{2) by g — (Rq,x) where Rq E SO(2) is the 



g = sin (9 cos 9 y J E SE{2). 

V oil 



cos 9 — sin 9 x 
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Center of mass Center of mass 




(a) Body frame is aligned with axes (b) Arbitrary position and orienta- 

of symmetry of the body tion of the body frame. 

Figure 1: Two different choices of the body frame for an elliptical two-dimensional rigid body. In 
both cases the origin of the body frame does not coincide with the center of mass. 



rotation matrix determined by the angle 9. Let (Vl, v 2 ) GK 2 be the linear velocity of the origin of the 
body frame written in the body coordinates, and denote by uj — 9 the body's angular velocity. They 
define the element £ in the Lie algebra se{2) given by 

/ -00 V! \ 

t = g- l g= [ uj v 2 \ ese(2). 

\ / 

Explicitly we have 

9 = uj, v\ — x cos 9 -\- y sin 9, v 2 = — xsin# + y cos 9. (2.1) 

For convenience, we will sometimes denote £ E se(2) as the column vector (cj,vi,V2) T £ M 3 . The Lie 
algebra commutator takes the form 

[ (cj, v 1 ,v 2 ) , (cj', v[, v' 2 ) ] 5e(2) = ( , v 2 J - uv 2 , uv[ - vxJ). 

The kinetic energy of the body is given by 

T# = - ((2]g + ma 2 + mb 2 )uj 2 + m{y\ + v 2 ) — mbuovi + macjv 2 ) , (2.2) 

where m is the mass of the body, (a, b) are body coordinates of the center of mass (see Figure [I]), and 
Is is the moment of inertia of the body about the center of mass. It is a positive definite quadratic 
form on se{2) whose matrix is the body inertia tensor 

(Tb + m(a 2 + b 2 ) —mb ma \ 
-mb m . (2.3) 

ma rn I 
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The fluid flow at a given instant. Consider now the motion of the fluid that surrounds the body. 
Suppose that at a given instant the body occupies a region BcK 2 . The flow is assumed to take place 
in the connected unbounded region U := R 2 \ B that is not occupied by the body. 

We assume that the flow is potential so the Eulerian velocity of the fluid u can be written as u = V<3> 
for a fluid potential <E> : U —¥ R. Incompressibility of the fluid implies that $ is harmonic, 

V 2 $ = on U. 

The boundary conditions for $ come from the following considerations. On the one hand it is 
assumed that, up to a purely circulatory flow around the body, the motion of the fluid is solely due to 
the motion of the body. This assumption requires the fluid velocity V$ to vanish at infinity. Secondly, 
to avoid cavitation or penetration of the fluid into the body, we require the normal component of the 
fluid velocity at a material point p on the boundary of B to agree with the normal component of the 
velocity of p. Suppose that the vector (X, Y) E R 2 gives body coordinates for p. The latter boundary 
condition is expressed as 

— = (yi - uY)ni + (v2 + cjX)n 2 , 
dn P edB 

where n = (^1,^2) is the outward unit normal vector to B at p written in body coordinates. These 
conditions determine the flow of the fluid up to a purely circulatory flow around the body that would 
persist if the body is brought to rest. The latter is specified by the value of the circulation k around 
the body as we now discuss. 

The potential $ that satisfies the above boundary value problem can be written in terms of the 
body's velocities v\,vi,w, in Kirchhoff form: 

$ = Vl fa + v 2 (j)2 + u\ + (f>o, (2.4) 

where fa, i — 0, 1,2, and x are harmonic functions on U whose gradients vanish at infinity and 
satisfy: 

dfa 
dn 



-n- i-1 2 ^ 

dB dn 



= Xn 2 - Yn 1: —— 

dB 9n 



= 0. 

dB 

The potential </>o is multi- valued and defines the circulatory flow around the body. The circulation k 
of the fluid around the body satisfies 



cp u 

JdB 



• dl = (p V0 O • dl, 

dB 



and remains constant during the motion. Figure [2] shows the streamline pattern of the flow determined 
by the motion of an elliptical body for different values of k. 

Disregarding the circulatory motion, the kinetic energy of the fluid is given by 

?> = f / \\V($-<f> )\\ 2 dA, 

where dA is the area element in R 2 and p is the (constant) fluid density. We have subtracted the 
circulatory part from the velocity potential, as it would give rise to an infinite contribution to the fluid 
kinetic energy [25] . 
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(a) k = 



(b) k = 2tt 



(c) k = 10 



Figure 2: Stream line pattern for an ellipse moving on the plane for different values of the 
circulation k. The major and minor semi-axes of the ellipse are A = 2,B = 1. The body frame is 
aligned with the principal axes of the ellipse and the velocity of the body satisfies v\ = v<i = 1 and 
uj = 5. 



By substituting (2.4) into the above, one can express Tjr as the quadratic form 

where M l jr, ICjr, i,j = 1,2, and Xjr are certain constants that only depend on the body shape. Explicitly 
one has (see [T7] for details), 

M% = -p[ ^dl = -p[ ^dl,i,j = l,2, X^ = -p/ x ^cfl 

Vjr = -p[ <f>i^dl = -p [ x^dl, i = 1,2. 
JdB on J dB dn 

These constants are referred to as added masses and are conveniently written in 3 x 3 matrix form to 
define the (symmetric) added inertia tensor: 



It ■= 

that defines Tjr as a quadratic form on se(2). 



1>T K-T 



Example. For an elliptic rigid body with semi-axes of length A > B > 0, the added masses and 
moments of inertia take on a particularly convenient form. The kinetic energy of the fluid is given by 

(see El) 
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where we have ignored the circulatory motion around the body. The corresponding added inertia tensor 
is thus given by 

l^ft \ 

I j 7 — pTr I B 2 ■ (2.6) 
V A 2 J 

We emphasize that this particular form of the added inertia tensor was derived under the assumption 
that the body frame is aligned with the symmetry axes of the ellipse [y — 0, r = s = in Figure [3] 
ahead). When this is not the case, the added mass tensor is more complicated, and in particular need 



not be diagonal, as is shown in (3.12). 



2.2 Rigid body motion in potential flow 

The total kinetic energy, T, of the solid-fluid system (excluding the circulatory motion) is the sum of 
the kinetic energy T# of the rigid body and the energy Tjr of the fluid. As both T# and Tjr are functions 
on TSE(2), so is the total energy T. In the absence of external forces or circulation, the Lagrangian 
C of the solid-fluid system is just the kinetic energy: C — T, and in this case, the motion of the rigid 
body describes a geodesic curve in SE{2) with respect to the Riemannian metric defined by C. 



In view of (2.2) and ( |2.5[ ), we can write the Lagrangian C — Tg+Tjr in terms of the linear and angular 
velocities of the body (written in the body frame) and this expression does not depend on the particular 
position and orientation of the body, i.e. is independent of the group element g = (i?#,x) E SE(2). 
Thus C is invariant under the lifted action of left multiplication on SE{2). This symmetry corresponds 
to invariance under translations and rotations of the space frame. The reduction of this symmetry 
defines Euler-Poincare equations on the Lie algebra se(2) or, in the Hamiltonian setting, the Lie- 
Poisson equations on the coalgebra se(2)*. The latter are precisely Kirchhoff's equations that are 
explicitly written below. 

The invariance of C allows us to define a function L on se(2), called the reduced Lagrangian. 
Explicitly, L is given by 

L(0 = \fK, 

where £ = (cj, vi, V2) T E se(2) = R 3 is a column vector and the matrix I is the sum of the inertia matrix 
Ib of the rigid body and the added masses and inertia Ijr of the fluid: I = I# + I jr. 

Since se(2) is isomorphic to M 3 and using the Euclidian inner product, we identify se(2)* ^ R 3 . A 
typical element fi is represented as a row vector fi = (fc,pi,p2)- The duality pairing between \i and an 
element £ = (a;, v\, V2) T of se(2) is given by 

(/i, f ) = (j£ = kcj + ptvx + p 2 V2- 
With this identification, the Legendre transform associated to L is defined as the mapping 

¥L : se{2) se(2)* 
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given by FL(£) = /i, where p — (I0 T - The components of p — {k,p\,P2) are explicitly given by 



fc = (1b + m(a 2 + 6 2 ) + Xjr)a; + (-raft + Kjr)vi + (ma + K,jr)v 2 , 
Pl = (-mb + K} T )u + (m + .M^i + X 12 ^2, (2.7) 
p 2 = (ma + JCjr)u + 7W 12 i;i + (m + TW^)^- 

In classical hydrodynamics k and (pi,P2) are known as "impulsive pair" and "impulsive force" respec- 
tively. 

The reduced Hamiltonian H : se(2)* — »> R is given by 

and the corresponding (minus) Lie-Poisson equations are p = adj_i^/x. Written out in component form, 
these equations are nothing but the Kirchhoff equations: 

k = v 2 p 1 -v 1 p 2 , , Q v 

Pl=^P2, P2 = 

where the velocities (cj, i?2) T and the impulses (fc,pi,P2) are related by the Legendre transformation 



(2.7). Finally, we remark that the motion of the body in space can be found from a solution of (2.8) 



by solving the reconstruction equations (2.1). 



2.3 Rigid body motion with circulation 

In the presence of circulation, the Kirchhoff equations on se(2)* have to be modified to include the 
Kutta-Zhukowski force mentioned in the introduction. This is a gyroscopic force, which is proportional 
to the circulation n. In this case, the equations of motion are referred to as the Chaply gin- Lamb 
equations [6} 117], and they are given by 

k = V2P1 - V1P2 - p(avi + /3v2), 
Pi = up2 - npV2 + pau, (2.9) 
P2 = —ojpi + npvi + p/3uj, 

The constants a and (3 are proportional to the circulation k and depend of the position and orientation 
of the body axes. They are explicitly given by: 



a = <b XV0o • rfl, 0= <b yV0 o • dl, (2.10) 

JOB JdB 

where, as before, (X, Y) are body coordinates for material points in dB. The Chaplygin-Lamb equations 
are discussed in detail in [3] and a derivation from first principles, using techniques from symplectic 
geometry and reduction theory, is given in [29J . 

One easily verifies that if the center of the body axes is displaced to the point with body coordinates 
(r, 5), so that the new body coordinates are X — X — r, Y — Y — s, then the circulation constants 
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relative to the new coordinate axes take the form a — a — (3 = (3 — sk. Thus, there is a unique 
choice of the body axes that makes these constants vanish. On the other hand, it is also desired to 
choose the body axes so that the total inertia tensor I is diagonal. For an asymmetric body, these two 
choices are in general incompatible, see e.g. |17j . 

For our purposes, the choice of body axes will be made to simplify the expression of the nonholo- 
nomic constraint. We therefore consider equations (2.9) in their full generality where a, f3 ^ 0, and I 
is not diagonal. This contrasts with the treatment in [29] where it is assumed that a — (3 = and with 
[3] where the complementary assumption, namely that I is diagonal, is made. 

It is shown in [12] that equations (2.9) are of Euler-Poincare type on a central extension of SE{2) 
and thus are Hamiltonian. 



3 The hydro dynamic planar Chaplygin sleigh with circulation 

We are now ready to consider the mechanical system of our interest which is the generalization of 
the hydrodynamic version of the Chaplygin sleigh treated in [llj to the case when there is circula- 
tion around the body. Recall that the classical Chaplygin sleigh problem (going back to 1911, [7]) 
describes the motion of a planar rigid body with a knife edge (a blade) sliding on a horizontal plane. 
The nonholonomic constraint forbids the motion in the direction perpendicular to the blade. In its 
hydrodynamic version, the body is surrounded by a potential fluid and the nonholonomic constraint 
models the effect of a very effective fin or keel, see [IT] . 

With the notation from section [2j we let {Ei, E2} be a body frame located at the contact point of 
the sleigh and the plane, and so that the Ei-axis is aligned with the blade (see Figure [5]). The resulting 
nonholonomic constraint is given by V2 = 0, and is clearly left invariant under the action of £i?(2), as 
it is solely written in terms of the velocity of the body as seen in the body frame. 

In the absence of constraints, the motion of the body is described by the Chaplygin-Lamb equations 



(2.9). In agreement with the Lagrange-D'Alembert principle, which states that the constraint forces 
perform no work during the motion, the equations of motion for the constrained system are 

k = V2P1 - V1P2 - p(av\ + /3^ 2 ), 

Pi = - KpV2 + POLUJ, (3-11) 

P2 = —ujpi + npvi + pjiuj + A, 

where the multiplier A is determined from the condition ^2 = 0. These equations have been shown to 
be of Euler-Poincare-Suslov type on the dual Lie algebra of a central extension of SE{2) in [12J. 

The total inertia tensor and the circulation constants a, (3. The behavior of the solutions of 



(3.11) depends crucially on the total inertia tensor I = Ig + Ijr that relates (^,^1,^2) to (cj, ^1,^2)5 and 



on the value of the circulation constants n and a, (3. 

The expression for Ig with respect to the body frame {Ei,E2} is given by (2.3) where m is the 
mass of the body, (a, b) are body coordinates of the center of mass, and I is the moment of inertia of 
the body about the center of mass. While this simple expression is independent of the body shape, 
an explicit expression for the tensor of adjoint masses Ijr can be given explicitly only for rather simple 
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geometries. A simple yet interesting one is that for an elliptical uniform planar body with the semi-axes 
of length A > B > 0. Assume that the origin has coordinates (r, s) with respect to the frame that is 
aligned with the principal axes of the ellipse and that the coordinate axes Ei E2 are not aligned with 
the axes of the ellipse, forming an angle v (measured counter-clockwise), as illustrated in Figure [3j 



Center of mass 




Contact point 



Figure 3: The elliptical sleigh. The blade makes an angle v with the major axis of the ellipse 
and the contact point has coordinates (r, s) with respect to the frame that is determined by the 
principal axes of the ellipse (in the diagram both r and s are negative). 



For this geometry, starting from the formula (2.6) for the added inertia tensor given in section 2.1 
one can show that 



11 jr pir 



^ — * } + s 2 (B 2 cos 2 v + A 2 sin 2 v) 
+ r 2 (A 2 cos 2 v + B 2 sin 2 v) - rs(A 2 - B 2 ) sin(2^) 

s(B 2 cos 2 v + A 2 sin 2 v) 
-\r(A 2 - B 2 ) sin(2^) 

— r(A 2 cos 2 v + B 2 sin 2 v) 
+ ±s(A 2 - B 2 ) sin(2^) 



s(B 2 cos 2 v + A sin 2 v) 
-\r(A 2 - B 2 ) sin(2^) 

B 2 cos 2 v + A 2 sin 2 v 
a2 ~ b2 sin(2^) 



-r(A 2 cos 2 v + B 2 sin 2 v) \ 
+ ^s(A 2 - B 2 ) sin(2^) 



sin(2v) 



A cos z v + B z sin 2 v 



. (3.12) 



The total inertia tensor, I = Ig + Ijr, of the fluid-body system is then given by 

/ I + m(a 2 + b 2 ) 2222 2222 

(A 2 — B 2 ) 2 000 99 — mb + pir[s(B cos v + A sin v) ma + prc[— r(A cos v + B sin v) 

+ ^ [ 2 4 + « (* cos v + A sin „) -l r (A 2 -B 2 )sin(2v)] + ± .( A 2 - S 2 ) sin(2.)] 

+r (A cos i; + B sin v) — rs(A — B ) sin(2-u)] 2 2 



— mb + p7r[s(B 2 cos 2 -U + A 2 sin 2 t>) 
-|r(A 2 - B 2 ) sin(2w)] 

ma + p7r[ — r(A 2 cos 2 i; + £> 2 sin 2 v) 
+ |s(A 2 - S 2 ) sin(2u)] 



m + p7r[S cos i; + A sin u] 



Ptt [ A2 ~ b2 sin(2t,)] 



^— — 2 B sin(2i;)j m + p-7r(A 2 cos 2 v + £? 2 sin 2 v) 



Notice that in the presence of the fluid, if v ^ n|, n £ Z, the coefficient I23 = I32 is non-zero. This 
can never be the case if the sleigh is moving in vacuum as one can see from the expression given for 



Ig in (2.3). The appearance of this non-zero term leads to interesting dynamics that, in the absence of 



circulation were studied in [IT]. For the above geometry, the circulation constants a and (3 defined by 



(2.10) are computed to be: 

a = K,(rcosv + ssint;), f3 — hz(—r sin v + scosv). 

Notice that, in the presence of circulation, the two constants, a and /3, can vanish simultaneously only 
if r = s = 0, that is, only if the contact point is at the center of the ellipse. 
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In the sequel we assume that the shape of the sleigh is arbitrary convex and that its center of mass 
does not necessarily coincide with the origin, which leads to the general total inertia tensor 



/ J -L 2 Lt \ 

I = \ -L 2 M Z , 

\ Li Z N ) 

and with arbitrary circulation constants a, /3. 



Detailed equations of motion. The constraint written in terms of momenta is v 2 — I 1 (k, p\ , p 2 
0. Differentiating it and using (3.11), we find the multiplier 



\T - 



A = — 



(i- 1 : 



33 



V 2 pi - VlP2 - p{0LV\ + (3v 2 ) 
Ujp 2 — KpV 2 + paUJ 
—Upi + KpVi + pf3(jj) 



where 



det(I) 



MN - Z 2 ZL X + NL 2 -ZL 2 - ML\ 
ZLi + NL 2 JN - L\ -L X L 2 - JZ 
-ZL 2 - MLi -L X L 2 - JZ 



JM-L\ 



A long but straightforward calculation shows that, by expressing uj,v\ and v 2 in terms of k,pi,p2. 



substituting into (3.11), and enforcing the constraint v 2 = 0, one obtains: 

Co = (L\uj + Zv\ + pa) (L 2 uj — Mv\) , 
vi = (Liuj + Zv\ + pa) (Juj - L 2 vi) , 



(3.13) 



where we set D = det(I)(I- 1 ) 33 = MJ - L\. Note that D > since I and I 1 are positive definite. 
Note as well that if a = we recover the system with zero circulation treated in [llj so from now on 
we assume 

The full motion of the sleigh on the plane is determined by the reconstruction equations which, in 
our case with ^2 = 0, reduce to 



UJ, 



X — Vi cos#, 



y = vi sin ( 



The reduced energy integral has 

H = ^ ( Juj 2 + Mv\ - 2L 2 ujv ± ) , 

and its level sets are ellipses on the (uj i?i)-plane. 

As seen from the equations, the straight line t — {L\uj + Zv\ + pa — 0} consists of equilibrium 
points for the system. Each of these equilibria corresponds to a uniform circular motion on the plane 
along a circumference of radius — . 
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Notice that if Z = L\ = the line £ of equilibria disappears. In fact, it is shown in [12J that 



equations (3.13) possess an invariant measure only for this choice of the parameters. In this particular 
case we obtain simple harmonic motion on the reduced plane uj,v\. 

For the sequel we will assume that Z and L\ are not both zero. We shall see that initial conditions 
with high energy yield asymptotic dynamics which is in agreement with our statement that there is no 
invariant measure in this case. 

A level set of the energy will intersect once, twice or never the line of equilibra £ depending on the 
value of H. One can show that there are two intersections if H > ho, only one intersection if H = ho, 
and zero intersections if H < ho, where 



with E = JZ 2 + 2ZL ± L 2 + ML\ > 0. 



Hence, the trajectories of (3.13) are contained in ellipses and they are of three types: 

1. For small values of the energy, < H < ho, we have periodic motion on the ellipses. 

2. For H = ho we have a homoclinic connection that separates periodic from asymptotic trajectories. 

3. For H > ho we have heteroclinic connections between the asymptotically unstable and stable 
equilibria on £. 

We introduce two other special energy values h\,h 2 > ho, for which the corresponding energy 
contour line intersects the equilibra line £ at the axes uj — and v\ = respectively. Namely, 

1 {pafM 1 (pa? J 

kl ~ 2 h2 ~ 2 "LP" ( } 

The phase portrait is illustrated in Figure [4] (a). 



Remark 1. In fact, the reduced system (3.13) can be checked to be Hamiltonian with respect to the 



following Poisson bracket of functions of uj,v\ 



{Fi, F 2 } :=-— {Liu + Zvi + pa) 



D \ duo dv\ dvi doj 

The invariant symplectic leaves consist of the semi-planes separated by the equilibria line £ and the 
zero-dimensional leaves formed by the points on £. 



3.1 The motion of the sleigh on the plane: qualitative description 

As for the reduced dynamics, the qualitative motion of the sleigh on the plane depends crucially on 
the energy value. We distinguish 7 regions on the reduced phase space as illustrated in Figure [4] (b). 
The regions are separated by the homoclinic orbit corresponding to the energy contour H — ho, the 
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(a) Reduced phase portrait. The (b) The 7 regions in the reduced phase 

stable and unstable equilibra are space, 
represented by filled and empty 
dots, respectively. 



Figure 4: Reduced phase portrait and the 7 regions in phase space (assuming that L\ > 0, 
a, Z < 0, and hi < Z^)- 



two pairs of heteroclinic orbits corresponding to the energy contour lines H — hi and H = /12, and 
the line of equilibra t. It is thus clear that the regions are invariant by the flow. We shall see that the 
qualitative motion of the sleigh is different in each region. The number of regions can change for special 
parameter values that make any of ho, hi, hi coincide. We shall discuss the behavior in the interior of 
the regions in the generic case where they are all different and under the assumption that hi < hi- 



Motion in Region I. For subcritical values of the energy, < H < ho, the reduced dynamics on 
(vi,(jj) is periodic, but the motion on the sleigh is generally not periodic. Let T be the (minimal) 
period of the reduced dynamics that depends on the energy value H in a way that will be made precise 
in section |4j Under the shift t —> t + T the coordinates of the sleigh undergo the increments 



Ax = I vi(t) cos 6 (t)dt, Ay = [ vi(t) sin0(t) dt, A6 = [ 
Jo Jo Jo 

These quantities completely determine the type of motion of the sleigh. 



u(t) dt. 



(3.15) 



Theorem 3.1. Let T be the period of a periodic solution to the reduced system (3.13) with H < ho, 
and let Ax = (Ax, Ay) and AO be defined by (3.15). Then, the motion of the sleigh on the plane is 

1. T -periodic if^eZ and Ax = 0. 

2. Unbounded if^eZ and Ax ^ 0. 
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3. Contained in a circle or an annulus if ^ where the motion 

(a) is periodic of period qT if ^ is a rational number written in terms of integers p,q as ^ = | 
in irreducible form, 

(b) fills up the annulus densely if ^ is irrational. 

Moreover, the above types of the motion depend only on the energy value H. In particular, they are 
independent of the initial configuration of the sleigh. 

Proof. Assume without loss of generality that at t = the body and the space reference frames coincide. 
Then 9(0) = 0, and the element 

cos(A6>) - sin(A0) Ax 
sin(A#) cos(A#) Ay \ e SE(2) 
1 

describes the shift in both reference frames after one period. Then the elements 

cos 9(nT) -sin 9 (nT) x(nT) 
g n =\ sin9(nT) cos9(nT) y(nT) ) e SE(2), neZ 
1 

describing the subsequent positions of the sleigh in the space frame, are the right translations on SE(2): 

9n+l 9nS. 

Following theorems of kinematics, each generic transformation of that kind is a rigid rotation about a 
fixed point Q on the plane (x, y) by the angle A9. This happens when ^ ^ Z. Then the trajectory of 
the contact point lies inside of an annulus or a circle with the center Q. The trajectory is periodic iff 
^ is a rational number, otherwise it fills the domain densely. 

In the special case ^ E Z the transformation g n g n +i is a parallel translation (if Ax ^ 0), 
which is obviously unbounded, or is the identity (when Ax = 0), then the trajectory is periodic. This 
leads to the items of the theorem. 

Finally note that the components of S depend only on the energy value and not on the particular 
initial conditions. □ 

We mention that the matrix S in the proof is the (right) monodromy matrix of the periodic linear 
system defined by the reconstruction equations. The global description of the invariant manifolds in 
the unreduced space, that generically are two-tori, can be studied with techniques similar to those 
developed in [8]. 

The types of motion indicated in items 2 and 3 of the above theorem are illustrated in Figure [5| 
The generic behavior of the sleigh corresponds to item 3 (b) (Figure [5] (c)). 

Equation (4.19) in section [4] below gives the explicit dependence of A9 on the energy. 



Now notice that in the case h > ho, the equilibria points (cj, v{) on the line £ = {L1UJ + ZV1 + pa — 0} 
correspond to circular motion of the sleigh in the direction determined by the sign of uo (clockwise if 
cj > and anti-clockwise if ou < 0). The circles describing by the contact point have radius \vi/cj\. 
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(a) Trajectory of the 
contact point with AO = 
— 47r, and Ax, Ay / 
(unbounded motion). 



(b) Trajectory of the con- 
tact point with AO — — ^ L 
(periodic motion in an annu- 
lus). 



(c) Trajectory of the contact 
point with AO = -^f for 
a finite time (the trajectory 
fills the annulus densely). 



Figure 5: Trajectory of the sleigh on the plane in region I (H < ho) for different values of AO. 



The exceptional cases when uj = or v\ — respectively correspond to motion along a straight 
line or to a spinning motion of the sleigh about the contact point that remains fixed. These cases 
correspond respectively to the energy values h\ and /i2- The motion of the sleigh on the plane when 
the energy attains the critical value ho will be discussed in section [4j 

Motion in Region II. The sleigh evolves asymptotically from one circular motion to another. The 
limit circles have different radii. Both of them are traversed counterclockwise (at the limit points in the 
region one has uj > 0). Moreover, the trajectory of the point of contact has two cusps corresponding 
to the two intersections with the axes v\ — on the Figure [4j 

Motion in Region III. The same behavior as in Region II, however, this time the trajectory has 
no cusps (since v\ never vanishes). 

Motion in Region IV. The sleigh again evolves asymptotically from one circular motion to another. 
The limit circles have different radius and are are traversed in opposite directions. The trajectory of 
the contact point has two cusps corresponding to the two intersections with the axes v\ — on the 
plane (vi,uj). 

Motion in Region V. The same behavior as in Region IV, but the trajectory has no cusps since v\ 
never vanishes. 

Motion in Regions VI and VII. The same behavior as in Region IV, however the trajectories 
have exactly one cusp corresponding to the intersection of the axis v\ — 0. 

Thus we see that the qualitative behavior of the motion of the sleigh on the plane can be sensitive to 
initial conditions. This contrasts with the properties of other completely solvable classical nonholonomic 
systems like the Suslov problem or the hydrodynamic Chaplygin sleigh without circulation. For these 



15 



(a) Region II. (b) Region III. (c) Region IV. 

Figure 6: Trajectories of the sleigh for initial conditions in regions II, III, and IV. 




Figure 7 : The trajectories for initial conditions in regions V, VI, and VII. 



systems one can obtain formulas for the parameters that determine the long term behavior of the 
dynamics, and such formulas are independent of the energy (see fTTj [10]). 

At the physical level, the motion of the sleigh can be understood as the interplay of two effects, 
the circulation of the fluid around the body, and the inertia of the body. For low energies (H < ho) 
the dynamics are driven by the circulation and produce periodic motion in the reduced space. For 
high energies (H > ho) the inertia of the body takes over and the asymptotic dynamics of the body 
on the plane resemble the motion of the hydrodynamic Chaplygin sleigh in the absence of circulation 
considered in fTT] . 

We finally stress that all the above types of motion also appear in the classical Appel-Korteweg 
problem on the rolling of a non-vertical disc. Namely the trajectories of the contact point of the disc 
admit the same classification (see [T| \19\ 19]). 

4 Explicit solution of the reduced system in se(2)* 



We will only treat the case — when the system (3.13) takes the form 



uj = ^- (L\uj + Zvi + pa) (-Mvi) , 

1 (4.16) 
v\ = — {Liu + Zvi + pa) (Ju) , 
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with D — MJ. The energy takes the diagonal form H — \{Jlo 2 + Mvf). The general case can be 
reduced to this one via an invertible linear change of variables as the following proposition shows. 



Proposition 4.1. There exists a 2 x 2 rotation matrix P £ 50(2) such that if lo and v\ satisfy (3.13) 
then (lo, vi) T :— P(lo, v\) T satisfy (4.16) with coefficients L\, Z, M, J, where J, M > satisfy MJ — D, 
and (L 1 ,Z) T := P(L 1 ,Z) T . 



Proof Define the restricted inertia tensor I c by 



J 

-L 2 



-L 2 
M 



Since I c is symmetric and positive definite, there exists P £ SO(2) such that PI C P T is a diagonal 
matrix whose entries J, M > are the eigenvalues of I c and satisfy MJ — detl c = D. 



Equations (3.13) can be written in vector form as 



d 
dt 



1 

D 



LO 
VI 



Li 
Z 



+ pa 



-1 

1 



LO 
V\ 



where (•, •) denotes the usual inner product in R 2 . Defining (lo,vi) t :— P(lo,v\) t , the above equations 
are equivalent to 



d_ 
dt 



LO 
V\ 



1 

D 



LO 
V\ 



Li 
Z 



+ pa ) P 



The result now follows by noticing that 



-1 

1 



IcP 1 







PlcP 1 



-1 

1 



-1 

1 



IcP 1 



J 
M 



LO 
V\ 



where we have used that P 6 50(2) in the first equality. 



□ 



The solution of the system (4.16) can be obtained by parameterizing level set H — h of the energy 
function H = \{Jlo 2 + Mv\) as 



u>(t) 



2h ( 2f(t) 



j \f{ty + i 



vi(t) 



2h (f(t)< 



m \f{ty + i 



Substitution into (4.16) gives a separable equation for f(t) of the form 

f(t) = Af(t) 2 + 2Bf(t) + C, 



(4.17) 



where the coefficients A, £>, C, depend on the entries of the inertia tensor I, on pa, and also on the 
energy value h. The form of the solutions for the above equation depends on the sign of the discriminant 
A := B 2 — AC. A long but straightforward calculation shows that 

A-i(£)(k-Ao), 
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and hence sign(A) = sign(/i — ho). Therefore the form of the solution varies depending on whether 
the energy h is higher than, smaller than or equal to the critical energy level ho- We will now give the 
explicit formulae for u and v\ in these three cases. To write our formulas in a compact and unified 
manner we introduce the coefficients: 



7 := 



-sign(aZ) 



V2D 



\/\h - h \, K 1 :-- 



Vh 



ZVJ \y/h~i + Vh 



E 



Zyfj 




Vh + y/hi 



3 ■ = 



1 



ho + y/hi 



(4.18) 



where h\ is given by (3.14), and we assume that the constants Li, Z, are non-zero. Notice that 7, K\, 



and K2 depend on the energy value h. 



The solution for < h < ho- In this case the discriminant A < and f(t) = ^(j t&n^t) — B) is a 
solution of (4.17). After some algebra this gives the following solution of (4.16): 

K2 cos 2 (7t) + K\ sin(7t) cos(7t) 



u(t) 



J \(1 + K\) cos 2 (7t) + 2K 1 K 2 sin(7*) cos(7t) + K\ sin 2 (7t) , 
( ) _ , . [2h ( {K\ - 1) cos 2 (7t) + 2K X K 2 sin( 7 t) cosfrt) + K\ sin 2 ( 7 t) \ 
V!()- sign(a >\j M\(\ + Kl) cos 2 + 2K X K 2 sin(Tt) cos(7t) + K\ sin 2 (7t) / ' 

We immediately obtain the following formula for the period T of the solutions: 

T= - = — 

7 



00 as h — > ho from the left as expected. 



h 



Notice that T 

We can also obtain a closed expression for the increment AO — J Q T u(t) dt. For this matter we note 
that a primitive of Lo(t) is given by (see e.g. [T5]): 



9(t) 



2h ( e^t) + e 2 (t) + e 3 (t) 

J \AKlKl + {l + Kl - KlY 



with 



9 l {t) = K 2 {l + K 2 l +Kl) 1 t, 



02$) = + ^2 - 1) ln ((! + ^2) cos 2 ( 7 t) + 2KiK 2 sin( 7 t) cos( 7 t) + K{ sin 2 ( 7 t)), 

6» 3 (t) = -2KiK 2 arctan(Ki tan^i) + K 2 ). 

Therefore, taking into account the change of branch in the expression for 0s(t), we find that, under the 
above assumption L 2 — 0, 



AO = 0(tt/7) - 9(0) 



2vr f2h ( K 2 {1 + K\ + K\) + 2K 1 K 2 \ 
7 V J V 4K 2 K 2 2 + (1 + X 2 - K\Y ) ' 



(4.19) 
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(In the general case L 2 7^ due to Proposition 4.1, A9 will be a linear combination of (4.19)) and of 
the integral J Q T vi(t) dt.) 

In particular notice that A9 = O((ho — /i) -1 / 2 ) as h —> ho from the left. This implies that the 
behavior of the sleigh on the plane is extremely sensitive to the energy values that are slightly smaller 
than ho. Namely, we have 

Proposition 4.2. For any 8 > there exists a strictly increasing sequence of energy values {frj^HeN 
satisfying ho~S < h^ < ho, for which the motion of the contact point in the plane is periodic. Similarly, 
there exists an increasing sequence of energy values {Z&j^HeN (respectively, {h^}keN) satisfying 
ho — S < h^ < ho (respectively, ho — 5 < h$ < ho), such that the motion of the sleigh is unbounded 
(respectively, quasiperiodic) . 



The proof follows from the results of Theorem 3.1 and the fact that |A0| 00 as h —> ho from the 
left. 



The solution for h = ho* In this case the discriminant A = 0, and (4.17) has the solution 
f[t) = —\{\ + B). After some algebra this gives the following solution of (4.16): 



/ , \ 2ho K 3 t + K 2 t 

= 2\/ — 



vi(t) = sign(aZ) 



J (l + K 2 )t 2 + 2K 3 K 2 t + K 2 ' 

2h^ (K 2 - l)t 2 + 2K 3 K 2 t + K 2 



M (l + K 2 )t 2 + 2K 3 K 2 t + K 2 



2 ' 



where it is understood that h — ho in the expression for K 2 . 

To analyze the trajectory of the contact point of the sleigh, consider the particular solution u(t — 
to), vi(t — to) with to = — 7X77§. In order to keep the notation as simple as possible, we will also denote 
this solution by uj(t), v\(t). We have 



( . _ /2fep ( C 2 K 2 t 2 - K 3 dC 2 t - K 2 K 2 \ 

Mt) = aria^M ( C 'W + ***** " K & 



M V d(C 2 t 2 + K 2 ) 



where d = 1 + K 2 and C 2 = K 2 
Notice that 



hm — -signfaZh/ — — - := it. 

hence the limit motions of the sleigh are along circles of the same radius \R\. Moreover, the circles are 
traversed in the same direction as both uo(t) and v\(t) have the same limits as t — » zboc. 

Proposition 4.3. The motion of the contact point in the plane in the case h — ho is bounded. 
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Proof. It is sufficient to show that the distance between the centers of the limit circumferences is finite. 
The centers coincide with limit positions of the point C of the sleigh with body coordinates (0,i?). 
Indeed, the components of the velocity of this point are given by 

x c = V c (t) cos(0(t)), yc = V c (t) sin(0(t)), 

where 

V c (t) = vi(t) - Bw(t) = sign(aZ) 



2h K 3 (4Ki + Ci)t 



M K 2 (C%t* + K$y 

that goes to zero as t — >> ±00. Therefore, the vector that connects the two centers has the components 

/OO Z'OO 
V c (t) cos 9(t) dt, Ay c = / V c (t) sin 9(t) dt, 
-OO J —OO 

where 0(t) = 0i(t) + 9 2 {t) + 3 (t) with 



n 



[2KKzK 2 fdt\ 
e 3 (t) = -^ — ^parctan^— j, 

00 being an integration constant. Notice that 0i,03 and Vc are odd functions of t whereas 9 2 is even. 
Therefore, using basic trigonometric identities, we get: 

/oo poo 
V c sin((9i) sin(0 2 ) cos((9 3 ) dt - / Vc cos(6>i) sin(0 2 ) sin(0 3 ) dt. 
-00 J —00 



The above integrals can be shown to converge using integration by parts. For example, using that 
At 



sin(0i) = 4z(— Cs cos(0i)) for a certain constant C3, the first integral rereads 



/oo /»oo 
y c sin(0i)sin(0 2 )cos(0 3 )dt = -C 3 cos(0i)F c sin(0 2 ) cos^)^ + C 3 / cos(9i)k(t) dt, 
-OO J — OO 

where fc(t) = |(Vb sin(0 2 ) cos(0 3 ))- It is seen that the two terms on the right are finite. Namely, 
the limits of the boundary terms are zero since Vc —> as t — > ±00. Next, the integral is absolutely 
convergent since \k(t)\ can be dominated by a function that decays as 1/t 2 for large \t\. 

The proof that Ayc is finite is analogous. □ 



The solution for h > ho. In this case the discriminant A > and two solutions of ( |4.17 ) should be 
considered. These are fi(t) = ^(— 7tanh(7t) — B) and f 2 (t) — ^(— 7 coth(7t) — £>), and correspond 
to the two heteroclinic orbits that, together with the equilibra, make up the energy contour. 

Starting out with fi(t) we find the solution to ( 4.16[ ) 

l~2h ( K 2 cosh 2 (7t) - K\ sinh(7*) cosh(7t) 



v\(i) — sign(aZ) 



J V (1 + K'i) cosh 2 (7t) - 2K\K 2 sinh(7t) cosh(7t) + K\ sinh 2 (7t) , ^ ^ 

2h ([K\- 1) cosh 2 (7t) - 2K X K 2 sinh(7*) cosh(7t) + K\ sinh 2 (7t) x 



M V (1 + K$) cosh 2 ( 7 t) - 2K X K 2 sinh( 7 t) cosh( 7 t) + K\ sinh 2 ( 7 t) 
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that is contained in the semi-plane L\oj + Zv\ + pa > 0. 
Whereas considering /2(f) one finds the solution to ( |4.16| ) 

/2fr / /Y 2 sinh 2 (<yt) - K\ sinh( 7 t) cosh( 7 t) 



J V (1 + #2) sinh 2 ( 7 t) - 2Kiif 2 sinh( 7 t) cosh( 7 t) + K\ cosh 2 ( 7 t) 



, . /2fr / (K| - 1) sinh 2 ( 7 t) - 2/Yi/Y 2 sinh( 7 t) cosh( 7 t) + /Yj 2 cosh 2 ( 7 t) 

t>i(t) = sign(ais)W — — 



M V (1 + sinh 2 ( 7 t) - 2iYiiY 2 sinh( 7 t) cosh( 7 t) + K\ cosh 2 ( 7 t) / ' 

that is contained in the semi-plane L±co + Zv\ + pa < 0. 

For the sequel, we only consider the branch of the solution corresponding to the expressions for 
u)(t) and vi(t) given by (4.20). 



Integrating ui(t) and denoting k — 2\/2h/ J , we get 

6{t) = J udt = e l +6 2 + 6 z , 
Ki + K 2 



(2Y1 + 2Y2) 2 + 1 *' 

e 2 -- Kl kI-kI-i 



2 7 [(i^i + K 2 f + 1] • [(Ki - K 2 f + 1] 

• In ([(Ki - K 2 f + l]e 4 ^ + [2(/T 2 - /T 2 + l)]e 2 ^ + \{K X + /T 2 ) 2 + 1]) , 

63 = ^ 1+ K 2 ) 2 + 2 ^ 1 -^) 2 + l] arCtan m ~ + + ^k [Kl ~ Kl + ^ 

We will assume that sign(aZ) = 1 in the expression for v\. Note that the radii of the limit circles 
are different: 

Hm vi = 1 [J (K 2 T ^i) 2 ~ 1 
r± 2 V M X 2 =F Ki 

Then, to evaluate the distance between their centers, introduce the "floating radius" 

P ^ = e lt + e -7f ' S ° that / 5 ( ± °°) = r ±' 

and the floating point of the sleigh whose coordinates in the body fixed frame are (0, p(t)). The limit 
positions of this point coincide with the centers of the limit circumferences. The components of the 
velocity of the point are 

x p = V p {t) cos 9(t), ij p = V p (t)sm0(t), 

where 

, , N , x , x [2h 4K 2 1 

V p (t) = Vl (t)-p(t)u(t) - ■ ' 1 



M Kl - K\ [{Ki - K 2 ) 2 + l]e 2 ^ + 2[/T 2 - /T 2 - 1] + \{K X + K 2 ) 2 + l] e - 2 ^ ' 
Now introduce the new variable 

1 pr 1 -tf 2 ) 2 + l]e 2 * 



2Ki 
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and the constant 



Then 



k:=K$- K\ + 1 



e 2 ^ = 



Eh + JZ 2 (2^/hhi + hi) 
2K x z 



(4.21) 



{K l -K 2 ) 2 + Y 



dz = ^-\(K X - K 2 f + l]e 27 * dt = 2jzdt, 



(4.22) 



Z + T^Z + 



+ 1 



(Kt - K 2 ) 2 + 1 



Kt 4K 2 
and, therefore, 

K 1 + K 2 1 



4K 2 



([(Kt - K 2 ) 2 + l]e 4 ^ + 2(i^ 2 2 - K\ + l)e 2 ^ + (i^ + K 2 ) 2 + l) . 



1\ = K 



j 2 = K 



i 3 = K 



In 



2Ktz 



(K 1 + K 2 ) 2 + 12 1 \(Kt -K 2 ) 2 + lJ ' 
iff -K\-\ 



2 7 [(Jd + fT 2 ) 2 + 1] • [(#i - # 2 ) 2 + 1] 

2K X K 2 

7 [(Ki + K 2 ) 2 + 1] • [(Kt - K 2 ) 2 + 1] 

2h [(Kt - K 2 ) 2 + 1] 1 



M 



arctan yz + 

27* _ 



A 2 

#r 1 4K 2 

us 



In I z 2 + -j^-z H — + 1 ) (up to an additive constant), 



2#i 
2ft 



2Ktz 



The components of the vector of the distance between the centers of the limit circles are 

/oo poo 
V p cos 0(t)dt, Ay p = / V p sin 0(t)dt 
-oo J —oo 



Then, in view of (14. 221), 



/2/i 2Ki 
~MK\-K 2 

l2h K r 



3) I 



oo 2 2 i K3 z 1 1 1 

COS0 



z 2 + ^z + £k + l 



cos 0(t) dt 

— dz, 



/V 2 17 " w * t- 7^* -r 4^2 

and similarly for Ay^. 

Under the next substitution 2 = tani/, the above integrals yield 

<sj2h/M K\ ±- (sin^) ±mi (cos^) T ^ ai+2a2 ^ exp [±ia^ arctan (t&nu + 27^)) 



Ax c ± i^yc 



l{Kl - Kl) 



r f'1\ 

-a ±iai I 
I) ° Jo 



ao 



2Kt 



a 2 



(Kt-K 2 ) 2 + l 

nKt(Kl-K\-\) 

2((Kt + K 2 ) 2 + l)((Kt-K 2 ) 2 + iy 



at 



(l + 2^T sin(2n) + ^2 cos 2 u 
K(Kt + K 2 ) 



du, 



(4.23) 



a 3 



2 7 ((tf 1 + K 2 ) 2 + 1) ' 

2nK x K 2 



7 ((Kt + K 2 ) 2 + l)((Kt-K 2 ) 2 + l) 
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The latter integrals represent a generalization of the classical Beta- function (see, e.g., [13J). Their 
product gives the square of the distance between the centers. 



Remark. We could not calculate the above integrals explicitly. Note, however, that in the limit 
h — >> oc, when the energy prevails over the circulation effect, due to (4.18) and (4.21), the integrals 
Q4.23D reduce to 

Axq ± iAyc = —Cq 1 ^ / (sin^) T ^(cos^) T ^ (cos(2ci u) =p isin(2ci u)) du. 

Z Jo 



c 



DZ 



ci 



DML\ 

~E ' 



co 



Note that the same reduction holds in case of zero circulation (ho — hi — in (4.21) implies k — 0), 
which was studied in [TT] . As was shown there, the product of the above integrals gives the square of 
the distance between the centers in the following closed form 



(Ax c ) 2 + (Ay c ) 2 



2tt£> 



c 



z2 \cj + c< 



cosh(^7r) — cos(ci7r} 
sinh(C7r) 



Conclusions and further work 



We presented one of the first examples of nonholonomic hydrodynamical system, which is related to the 
design of underwater vehicles. The nonholonomic constraint can be interpreted as a first approximation 
model for a fin. From the mathematical point of view, our example is remarkable since both asymptotic 
and periodic dynamics coexist in the reduced phase space. It has been observed that the value of the 
energy is a crucial parameter in the qualitative behavior of the body on the plane. To our knowledge, 
this type of phenomenon is quite rare in nonholonomic dynamics. (A similar behavior occurs in the 
classical Appel-Korteweg problem on the rolling disc [HE].) 

For the future, we intend to study the motion of the hydrodynamic Chaplygin sleigh coupled to 
point vortices in the fluid [21]. The equations of motion for interacting point vortices and rigid bodies 
(without nonholonomic constraints) were recently derived in [261 E] and since then there have been 
significant efforts towards discerning integrability and chaoticity [22} [24] and towards uncovering the 
underlying geometry of these models [28]. We plan to couple the nonholonomic Chaplygin sleigh with 
one or several point vortices in the flow, taking these models as our next starting point. 
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